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Abstract. Gravitational waves can act like gravitational lenses, affecting 
the observed positions, brightnesses, and redshifts of distant objects. Exact 
expressions for such effects are derived here in general relativity, allowing for 
arbitrarily-moving sources and observers in the presence of plane-symmetric 
gravitational waves. At least for freely falling sources and observers, it is 
shown that the commonly-used predictions of linear perturbation theory can be 
generically overshadowed by nonlinear effects; even for very weak gravitational 
waves, higher-order perturbative corrections involve secularly-growing terms 
which cannot necessarily be neglected when considering observations of sufficiently 
distant sources. Even on more moderate scales where linear effects remain at least 
marginally dominant, nonlinear corrections are qualitatively different from their 
linear counterparts. There is a sense in which they can, for example, mimic the 
existence of a third type of gravitational wave polarization. 


1. Introduction 

Some of the most important potential signatures of gravitational waves are associated 
with their effects on the propagation of light. Collections of null rays can be 
deflected, sheared, delayed, or otherwise altered as they travel through a gravitational 
wave. Indeed, most contemporary attempts to observe gravitational waves rely on 
measurements of the relative time delays which accumulate as light travels between 
material bodies. This is particularly clear for interferometric detectors [1], where 
one or more beams of light are circulated between collections of mirrors and then 
recombined to reveal their relative phases. Efforts to detect gravitational waves 
using pulsar timing arrays [2] exploit similar principles, but instead make use of 
time intervals observed on the earth between radio bursts emitted by distant pulsars. 
Besides temporal effects such as these, gravitational waves can also affect observations 
of an object’s sky location, brightness, shape, and so on [3, 4, 5, 6, 7, 8, 9]. 

Almost all prior discussions of these phenomena have been perturbative, involving 
calculations which are valid only through first order in the gravitational wave 
amplitude (see, however, [3, 10, 11, 12, 13]). This has been justified, at least 
implicitly, by the minuscule size of even these lowest-order terms: In most cases 
of astrophysical interest, the gravitational wave strain amplitude e is much smaller 
than unity. Enormous technological effort is required to detect such waves at all, and 
waveform measurements which are accurate to more than a handful of decimal places 
cannot be expected for quite some time. In this context, it might appear reasonable 
to dismiss higher-order corrections as uninterestingly-small. One of the goals of this 
paper is to show that such reasoning can be misleading. Even if a dimensionless 
observable associated with a gravitational wave of amplitude e ^ 1 is bounded by 
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Figure 1: Schematic of a physical system which could correspond to the model 
considered in this paper. A “source” S emits electromagnetic radiation which is 
viewed by an “observer” O. In between these objects, spacetime is assumed to be 
approximately flat except for a nearly-planar gravitational wave. This wave may be 
generated by a distant binary, although all considerations here are restricted to the 
boxed region, and are therefore indifferent to the precise nature of the wave generation 
mechanism. 


e X (number of order 1) in linear perturbation theory, higher-order corrections are not 
necessarily bounded by x (another number of order 1). The coefficient in front of 
the term can instead grow enormously with the distance between a light source and 
its observer, implying that nonlinearities may be significant even when considering 
observations of very weak gravitational waves. Nonlinear effects also tend to have 
very different observational signatures from their lower-order counterparts, further 
increasing their potential detectability. 

Although it does not appear to have been previously pointed out in this 
context, the existence of large higher-order corrections is well-known in many types of 
perturbative calculations. A simple example is provided by the Mathieu-type equation 

C(w) + ^eC(M) cosu = 0. (1) 

If u denotes an appropriate phase coordinate, C^(u) may be shown to describe a 
particular metric component associated with a linearly-polarized, “monochromatic” 
gravitational plane wave with strain amplitude e. Moreover, the coordinate system 
where this is true is constructed such that there is a sense in which electromagnetic 
observations of distant objects have properties which can be read off directly from 
^^(u). Solutions to (1) therefore serve as a convenient proxy for understanding 
nonlinear effects associated with monochromatic gravitational waves in general 
relativity. Assuming e <C 1 while adopting convenient initial conditions, 

^^(u) = 1-1- ecosu + ie^(3cos^ u — u^) + O(e^). (2) 

8 

The magnitude of the second-order term in this expansion clearly overtakes the first 
when |u| ^ 3> 1, signaling that the linear approximation fails for large |u|. This 

occurs no matter how small e may be; weaker amplitudes merely delay such problems 
to larger scales. 

We show that similar effects arise for a variety of gravitational wave observables in 
general relativity, thus implying that the results of linear perturbation theory cannot 
necessarily be applied on large scales. As a model, geometric optics is considered 
in the presence of a plane-symmetric gravitational wave. This may be viewed as an 
idealization of the system illustrated in Figure 1, where observations are performed 
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sufficiently far from a gravitational wave source that the curvature of the wavefronts 
may be neglected. Similar models are common (though restricted only to first-order 
metric perturbations) in, e.g., descriptions for how gravitational waves can affect 
pufsar timing measurements [2, 14]. 

Despite this, real astrophysical observations cannot rely solely upon plane wave 
calculations. Deviations from planarity, waves propagating in multiple directions, 
non-radiative metric perturbations, and other effects could all be signihcant in 
observationally-relevant systems. Although calculations which take into account 
many such effects have been performed through first post-Minkowskian order [7, 15], 
the optical characteristics of nonlinear gravitational waves have been essentially 
unexplored in this context. Moreover, it would likely be difficult to understand 
the implications of any such calculations even if they did exist; the known first- 
order expressions are already extremely complicated in their most general forms. 
Plane waves are, by contrast, sufficiently simple that their physical effects can be 
thoroughly explored even in the nonlinear regime. At the same time, these waves 
remain sufficiently complicated to be interesting, and also to capture much of the 
relevant physics. Results obtained using plane wave descriptions may therefore be 
useful in the construction of specific hypotheses whose generality can later be tested 
using more complicated models. The technical details of the plane wave problem can 
also be used to suggest potential simplihcations in more general calculations. 

A completely separate motivation for considering plane wave spacetimes follows 
from a mathematical device known as the Penrose limit [16, 17, 18]. This provides 
a sense by which the geometry near any null geodesic in any spacetime is equivalent 
to the geometry of an appropriate plane wave. It can be interpreted as a statement 
that the metric in a small region around any sufRciently-relativistic observer in any 
spacetime is equivalent to that of an “effective plane wave.” Although we make 
no attempt to prove it, the Penrose limit suggests that (at least some types of) 
observations performed by ultrarelativistic observers can in general be reduced to 
analogous observations in effective plane wave spacetimes. 

Section 2 reviews gravitational plane waves in general relativity, first from the 
viewpoint of perturbation theory, and then as exact solutions to Einstein’s equation. 
Relations between these two perspectives and their relative advantages are described 
in detail. Next, Section 3 considers the physical consequences of plane wave spacetimes 
by deriving exact time delays, frequency shifts, observed sky positions, area distances, 
and luminosity distances. With appropriate identifications, some of the resulting 
expressions are only slightly more complicated than their linearized counterparts. 
Formal perturbative expansions are nevertheless derived in Section 4 and then applied 
to specific examples in Section 5. 

Notation and conventions 

The metric signature here is -1-2, c = G = 1, the Riemann tensor satisfies Rabc‘^^d = 
2V[oVf,]a;c for any Wc, Latin letters a,b,... denote abstract indices, Greek letters 
/r, iz,... denote four-dimensional coordinate indices, and i,j,... are used as two- 
dimensional coordinate indices associated with directions transverse to the background 
gravitational wave. When convenient, transverse coordinate components are indicated 
using boldface symbols without indices [e.g., jij = 'Yijw’‘v^ = w'^'yv, and 

tr 7 = Sijjij]. The one exception where the boldface symbol doesn’t correspond 
to its component counterpart is the 2x2 identity matrix {I)ij = Sij. Lastly, overdots 
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are used to denote derivatives with respect to a phase coordinate u, so 7 ^^ = d'yij/Au. 


2. Gravitational plane waves 

Before describing how light propagates in a gravitational plane wave, it must first 
be explained precisely what a plane wave is. Although the concept is clear for a 
scalar field in flat spacetime, subtleties arise when considering curved geometries or 
fields with nontrivial tensorial structure. This difficulty is reflected in part by the two 
distinct perspectives on gravitational plane waves—one perturbative and one exact— 
which are common in the literature. While both of these perspectives are individually 
well-known, the relations between them are not. This section clarifies the situation, 
and also remarks on special types of plane waves. 


2.1. Approximate plane waves 


Almost every textbook on general relativity discusses gravitational waves as linear 
perturbations on a flat background spacetime [19, 20]. To review, suppose that there 
exist coordinates (t, x^,x^, z) in which the metric components can be approximated by 
5/ji/ = ? 7 /i!y + e^/ji/ + 0(e^), where e ^ 1 is a dimensionless expansion parameter and 77 ^ 1 , 
is the Minkowski metric. It is also typical to adopt the transverse-traceless, or “TT” 
gauge, in which case hfj,t = = 0. These constraints on the metric 

perturbation can always be imposed in connected regions of spacetime where the 
linearized vacuum Einstein equation holds [21]. Imposing TT gauge in such a region, 
the first-order vacuum Einstein equation reduces there to the ordinary flat-spacetime 
wave equation for each coordinate component of the perturbation: = 0 . 

The solutions of interest here represent plane-symmetric gravitational waves 
propagating in vacuum. If the spatial projection of a plane wave’s direction of 
propagation is identified with djdz, the only components of the metric perturbation 
which might not vanish are hij, where i = 1,2. This 2x2 symmetric matrix must 
be trace-free, and can depend only on the “phase coordinate” u = {t — z)l'</2. The 
first-order line element for an arbitrary plane wave propagating in the z-direction is 
therefore 


ds^ = —dt^ -I- [5ij + ehij{u)\dx^dx^ -I- dz^ -I- O(e^), 


( 3 ) 


where 


h{u) 


h+{u) h^{u) \ 

hx{u) -h+{u) J’ 


( 4 ) 


and hx(u) represent the waveforms associated with the 

polarization states. A plane wave is said to be linearly polarized if hx 
made to vanish via some constant rotation of the x^, x^ coordinates. 


‘W” and “x” 
(or /i+) can be 


2.2. Exact plane waves 

The exact theory of gravitational plane waves is often presented very differently from 
its linearized counterpart. One way to motivate an exact plane wave solution in 
general relativity is to first search for those geometries which share—independently of 
Einstein’s equation—the symmetries associated with more familiar plane waves in flat 
spacetime. Consider, for example, a scalar field with the form / = f(t — z) in special 
relativity. This is clearly symmetric with respect to the two spacelike translations 
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djdx^ and the single null translation djdt + djdz. Less obviously, scalar plane waves 
are also preserved by 


V dz ) 


+ 




( 5 ) 


two Killing vector fields which generate rotations in the planes combined with 
boosts along the a:* directions. It follows that scalar plane waves in flat spacetime are 
preserved by at least five Killing vector fields associated with the geometry through 
which they propagate. The same symmetries also preserve electromagnetic plane 
waves in flat spacetime, thus motivating a gravitational plane wave in four dimensions 
as a curved spacetime which admits at least five linearly-independent Killing vector 
fields [22]. The resulting metrics are most commonly-stated in terms of the so-called 
Brinkmann coordinates {U,V, , X'^), in which case 

ds^ = -2dUdV + Hij{U)X^X^dU^ + (dXi)^ -p {dX'^f. (6) 


Using this as an ansatz, the exact vacuum Einstein equation reduces to 


trH = 0. 


( 7 ) 


Any 2x2 symmetric trace-free matrix Hij{U) therefore describes an exact plane wave 
in vacuum general relativity. We call this matrix the Brinkmann waveform. It has 
a direct geometric significance in the sense that the only independent, non-vanishing 
components of the Riemann tensor are 


RjJiUj — y (^) 

where i,j refer to coordinate components associated with the two X* coordinates. 
Furthermore, the two independent components of Hij correspond to the two possible 
polarization states of a gravitational wave in vacuum general relativity. Up to 
coordinate ambiguities in the construction of RuiUjy it follows from (8) that 
the Brinkmann waveform may be obtained using only local measurements. Such 
ambiguities are minor, and can be taken into account using only three constants Cb, 
Cr, and Ct [23]. Explicitly, two waveforms are physically identical if and only if they 
can be related via the complex replacement rule 

iLn(C/) + fiLi 2 (U) ^ cl[H^MU - q) + iH^ 2 {cbU - cO]U'=U (9) 

where C{, ^ 0 describes a constant boost along the direction of propagation, c* a 
constant translation of the phase coordinate U, and Cr a constant rotation. This 
is sufficiently simple that it is typically evident by inspection whether or not two 
Brinkmann waveforms describe the same physical system. 

Physically, a plane wave described by the metric (6) propagates in the null 
direction £“ = djdV. Up to an overall constant, £“ is the unique nonvanishing vector 
held which is both null and covariantly constant in any curved region of a plane wave 
spacetime. Noting that £a = —'^aU, these constraints define the U coordinate up to an 
overall affine transformation. It may be interpreted as the phase of the gravitational 
wave. Interpretations of the remaining Brinkmann coordinates follow by noting that 
they form a kind of Fermi normal coordinate system whose “origin” is the null geodesic 
V = X^ = 0 [17]. 

Plane waves are extremely simple when described in terms of Brinkmann 
waveforms. It is clear from (7) that if Hij{U) and Hl^iU) are any two vacuum 
waveforms, HijiU) -I- H[^{U) is also a vacuum waveform. This provides a sense by 
which plane waves in general relativity satisfy exact linear superposition; there is no 
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nonlinearity to speak of. Superpositions of this type are special cases of the more 
general result that the vacuum Einstein equation is linear for all metrics within the 
Kerr-Schild class. More precisely, suppose that there exists some null la and some 
scalar % such that pab = Pab + 'Hlalb is an exact solution to the vacuum Einstein 
equation for some vacuum “background” pab- If p'ab — Pab + H' I alb is a second exact 
solution, the metric 5 "^ = pab + (H + 'H')lalb must be an exact solution as well [24]. 
For the plane wave case of interest here, the Kerr-Schild decomposition is recovered 
by letting la = la = —^aU and H = , in which case pab is hat. Brinkmann 

waveforms are very special, and the linearity of Einstein’s equation which they make 
manifest is not at all apparent if plane waves are parametrized using different variables. 
Nevertheless, the optical observables discussed below are more conveniently described 
in terms of different waveforms which i) generalize the perturbative hij appearing in 
(3), and ii) do not satisfy linear superposition. 

2.3. Relating the exact and approximate descriptions 

Although the TT-gauge metric (3) describes, at least approximately, the same physical 
system as the Brinkmann metric ( 6 ), the former expression is not a trivial linearization 
of the latter. Understanding how these two descriptions relate to one another requires 
an appropriate coordinate transformation. The TT gauge makes sense only at first 
perturbative order, so a non-perturbative generalization of this gauge must be sought. 
We choose to employ those transformations which preserve, in an appropriate sense, 
the TT-gauge property that objects at “fixed spatial coordinates” move on timelike 
geodesics. Additionally, we require that the planar symmetry of the metric be manifest 
in the sense that there exist two coordinates cc®, interpreted as parametrizing the 2 - 
surfaces transverse to the gravitational wave, such that the vector fields djdx'" are both 
spacelike and Killing. These constraints are chosen not merely to facilitate a simple 
translation between the perturbative and non-perturbative viewpoints. Much more 
importantly, it is shown in Section 3 below that the “TT-like” coordinates which they 
define are particularly well-adapted to describing non-perturbative optical observables 
in the presence of an arbitrarily-strong gravitational wave. 

The appropriately-generalized TT coordinates are denoted here by (u,u, 
and are known in the literature as Rosen coordinates. They are related to their 
Brinkmann counterparts {U,V, X^, X'^) via 

u = u, x = r\U)X, v = V-^[iiU)rHU)],,X^X^, (10) 

where ^ij(U) is any nonsingular 2 x 2 (not necessarily symmetric) matrix satisfying 

iiU)=H{U)m, (€^Ofol= 0 . ( 11 ) 

Any with these properties is essentially a Jacobi propagator: Contracting it on 
the right with an arbitrary constant vector results in the nontrivial components of 
a solution! to the geodesic deviation (or Jacobi) equation [18]. Given any such 
propagator with the specified properties, the exact plane wave metric ( 6 ) becomes 

ds^ = —2dudv + ^ij{u)d.x''Ax\ ( 12 ) 

X The geodesic deviation equation may be written as ijj = —Rip in any spacetime, where -ip denotes 
the four-dimensional deviation vector resolved into parallel-propagated tetrad components, H is a 4x4 
matrix of similarly-decomposed Riemann components, and overdots represent ordinary derivatives 
with respect to an afhne parameter [25]. 
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= 7fo)(w) = ^k^{u)^kJ{u). 


(13) 


The implied summation in this last equation is understood to be trivial in the sense 
that or equivalently 7 = in matrix notation. Regardless, the 

Rosen line element (12) naturally splits into a longitudinal component —2dudu and a 
transverse component 'yij(u)dx^dx^. The longitudinal portion is flat and Lorentzian, 
while the transverse portion is associated with the Riemannian 2-metric ^ij(u). We 
refer to this 2-metric as the Rosen waveform, and note that it depends only on the 
phase u of the gravitational wave. It is related to the Brinkmann waveform Hij{U) 
via (11) and (13). 

That the Rosen coordinates are, as claimed, a type of geodesic normal coordinate 
system may be verified by noting that every fixed-x* 2 -surface contains a timelike 
geodesic. Moreover, each such surface actually contains a 1-parameter family of 
timelike geodesics related by longitudinal boosts. It also contains a null geodesic 
along which v = (constant). These statements can be made more intuitive in terms 
of the quasi-Cartesian coordinates 






(14) 


defined by analogy with standard null —>■ inertial transformations in flat spacetime. 
The vector d/dz is now a spatial projection of the wave propagation direction 
Additionally, any worldline which remains at fixed {x'',z) is a timelike geodesic. 
More generally, all curves satisfying dx^/dt = 0 and |dz/dt| = (constant) < 1 are 
timelike geodesics. Worldlines constrained by dx^/dt = 0, |d 2 ;/dt| = 1 are instead 
null geodesics. These statements imply that a large class of geodesics look trivial in 
Rosen coordinates. Geodesics which are not in this class can appear quite complicated, 
however. A generic timelike geodesic may be specified by choosing initial conditions 
for its transverse components at some fiducial phase u = Uo, as well as a constant 
A > 0 which describes motion longitudinal to the wave. Defining the phase average of 
the inverse 2 -metric by 


(7 = 


1 

U — Uo 



7 ^('r)dT, 


(15) 


an arbitrary timelike geodesic may be shown to have the coordinate parametrization 


x(u) = x(uo) + (u - uo)('y ^}'y(uo)x(uo), (16) 

v(u) = v(uo) + (u - Uo) + ^['r(uo)('r~^)'r(uo)]ijX^(uo)x^(uo)] . (17) 


Geodesics which remain at fixed (x\ z) satisfy i*(uo) = 0 and A = 1. 

That the Rosen metric directly generalizes the TT-gauge plane wave (3) follows 
from applying (14) to (12), which produces the exact line element ds^ = —dt^ -I- 
'Yij(u)dx‘dxl _|_ Standard perturbative results are therefore recovered if 

lij = + O(c^) (18) 


and tr/i = 0. It is shown in Section 4.1 that such expansions do indeed arise when 
considering smooth 1-parameter families of plane wave spacetimes. While the vacuum 
Einstein equation is both linear and algebraic in terms of Hij [cf. (7)], it is a nonlinear 
differential equation when expressed in terms of 7 ^-. This means that the O(e^) terms 
neglected in (18) are generically nontrivial. They may be viewed as corrections due 
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to the higher-order Einstein equation, or alternatively as higher-order solutions to a 
family of geodesic equations. 

The second of these perspectives follows from the interpretation of the spatial 
Rosen coordinates as a lattice of timelike geodesics. The relative displacements of 
these geodesics are encoded in Remarkably, this matrix also encodes in a simple 
way many properties of the null geodesics which are so central to the calculations 
of geometric optics. Plane wave spacetimes are sufficiently simple that most of the 
information required to characterize the propagation of light can i) be embedded in 
the Rosen-coordinate metric components, and ii) this coordinate choice is one for 
which considerable intuition and experience already exists (at least perturbatively). 
For these reasons, the remainder of this paper works almost exclusively in Rosen 
coordinates. 

More precisely, it is assumed that a particular Rosen coordinate system has been 
hxed. While the Brinkmann waveform Hij describing a particular plane wave is unique 
up to the relatively simple transformations (9) , there exist many physically-equivalent 
Rosen waveforms ^ij. Indeed, it is clear from (10) that relations between these two 
types of waveforms are nonunique and also nonlocal; different solutions to ( 11 ) are 
possible, and different resulting solutions for generically define different ^ij via 
(13). This freedom corresponds to the ability to choose different collections of timelike 
geodesics as coordinate markers. It can be described more precisely as the set of all 
solutions to for which is symmetric (a constraint which is true everywhere 

if it is true anywhere), and for which 7 = is nonsingular, modulo those solutions 
which preserve Fixing a particular Hij while choosing a fiducial phase uq, any 
two such solutions to ( 11 ) must be related by 


I'(m) = |(u)|| ^(mo)4'(mo) + (m - uo)(7 V{uo)i'{ uq) - (uq)^'{uq) |. (19) 


The complexity of this relation makes it clear that if one is presented with two 
Rosen waveforms 7 ^^-, 7 b, it not necessarily obvious by inspection whether or not 
they describe the same spacetime. It is nevertheless possible to check equivalence 
by computing some ^ij and which generate the appropriate waveforms, and then 
determining if and can be related by the transformation (9) [cf. (11)]. 

Although it is rarely discussed, the same ambiguities regarding representations of 
gravitational waves arise even in the linearized theory. As a possible counterpoint, 
one might object that vacuum, first-order, TT-gauge metric perturbations on 
a flat background are known to be invariant with respect to first-order gauge 
transformations. Such statements depend, however, on global assumptions such as 
asymptotic falloff (cf. Section 2.3 of [21]). If a spacetime is taken to be a literal 
plane wave, there is no such falloff. If, more realistically, a plane wave is used only 
as a model intended to be valid in a finite region of spacetime, the true asymptotic 
boundary conditions are external to the modeled region and therefore unusable. From 
either perspective, the familiar uniqueness of TT-gauge metric perturbations is lost. 
If (18) holds, the transformations hij{u) —>■ hij(u) + Cij + udij are easily shown to 
describe the same geometry through 0 (e) for any constant trace-free matrices Cy , dij. 

It is assumed here that a particular Rosen coordinate system has been chosen such 
that det^ > 0 everywhere of interest. This ensures that a single coordinate patch can 
be used in all calculations, and also that the orientations of the a;* and X* coordinates 
are identical. As with any geodesic-type coordinate system, gravitational focusing 
generically causes Rosen coordinate systems to break down on sufficiently large scales. 
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Restricting to a single coordinate patch therefore implies an upper bound on the 
maximum distances over which the calculations described here can be applied. This 
bound is, for example, of order in the case of a monochromatic gravitational 

wave with angular frequency uj and strain amplitude e. Optical effects involving larger 
scales are discussed in [3]. 

2.4- Rosen metrics: Non-perturbative considerations 

Suppose that a particular vacuum plane wave has been fixed by prescribing a trace- 
free matrix Hij{u). This is equivalent, via (8), to prescribing the wave’s curvature. 
Rosen waveforms then follow from (11) and (13), implying that ^ij{u) is essentially the 
square of a matrix which describes the displacements of coupled parametric oscillators 
attached to “springs” whose squared natural frequencies§ are proportional to Hij (it). 
Similar equations arise throughout physics and engineering, and a variety of methods 
have therefore been developed to understand them. 

A large class of physically-relevant waveforms can be classified as either 
“burstlike”—having a large amplitude only for a short time—or “continuous,” 
implying approximate periodicity. The curvature associated with a burst might be 
further idealized as vanishing completely whenever u lies outside a bounded connected 
region C C K. The spacetime can then be viewed as a curved region sandwiched 
between two flat connected regions J^±. Such spacetimes are often described as 
“sandwich waves.” In either of the flat regions where Hij = 0, the general solution to 
(11) is 

^{u) = ajr^+ubjr^, ( 20 ) 

The two matrices bjr^ must be constant and must satisfy i(ijrj_bj^j^)[ij] = 0 and 
det(ajr^ -b ubjr^) > 0 everywhere of interest. 

The trivial case of a spacetime which is everywhere flat is conventionally described 
by a metric in which = (5^-, implying that must be a constant orthogonal 
matrix. This is not required, however. Regions of flat spacetime can also be covered 
by coordinate systems whose “waveforms”—computed as squares of expressions 
like (20)— grow quadratically with u. Recalling (6) and (10), this corresponds to 
using non-comoving geodesics as coordinate markers. While such possibilities are 
an uninteresting complication in spacetimes which are globally flat, they can be 
unavoidable when considering bursts of gravitational waves. If coordinates are chosen 
such that "fij = Sij in one of a sandwich wave’s flat regions, the 2-metric in the 
other flat region cannot be prescribed arbitrarily, but must instead be computed by 
integrating (11) through the intervening burst of gravitational waves. Cases where 
resulting 2-metric is not equal to 6ij physically correspond to gravitational wave bursts 
which impart permanent displacements (and possibly kicks) to initially-comoving test 
particles. This is known as the gravitational memory effect, and is discussed further 
in Sections 4.1.1 and 5.2 below. 

Plane waves which cannot be modeled as bursts, but whose curvature is instead 
periodic in it, may be understood using different methods. In these cases, ^ 
corresponds to a set of coupled Hill-like equations. Although some methods exist for 
finding exact solutions to Hill equations [26, 27], these are rather limited. General 
properties of periodic waves may instead be understood using Floquet theory [28]. 

§ Einstein’s equation requires that the eigenvalues of Hij have opposite signs, so the instantaneous 
frequencies in this analogy are not necessarily real. 
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(a) Diagram illustrating the observa¬ 
tion of a pointlike source S by an ob¬ 
server O. The observation event po is 
connected to the emission event ps by 
a light ray whose tangent can be com¬ 
puted as the gradient of Synge’s func¬ 
tion a. Source and observer 4-velocities 
are denoted by . 


(b) A two-arm interferometer with an 
observer O flanked by two mirrors 
5 i, 2- The relative phases of the 
light compared at po determine the 
difference Sto in their emission times. 
This system may be viewed as a 
composite the one in part (a). 


It follows from this that even though and "fij = ^ki^kj do not typically share the 
periodicity of Hij, knowledge of ^ij in only one oscillation period can be used to extend 
it everywhere. Although we shall not do so here, such methods can be used to greatly 
enlarge the validity of perturbative calculations. Another interesting consequence of 
Floquet theory is that it can be used to find which waves are “unstable” in the sense 
that they admit exponentially-growing Such dramatic effects require very large 
gravitational wave amplitudes, however. 

3. Non-perturbative optical observables 

How, at least in principle, might a gravitational wave be detected? Perhaps the most 
obvious observable associated with any spacetime is its curvature, and it follows from 
(8) that knowledge of Rabc^ for any gravitational plane wave immediately implies 
knowledge of its Brinkmann-type waveform. Local experiments may therefore be used 
to directly determine Hij{u) for all phases u in which measurements are performed. 
Not all interesting observables are local, however. Pulsar timing techniques, for 
example, use radio bursts which are expected to travel through many gravitational 
wave cycles before reaching the earth. Observed properties of these bursts are therefore 
modified by a kind of integrated curvature which depends nontrivially on Hij. 

This section considers more generally those nonlocal observables which describe 
the appearance of a luminous source in the presence of an intervening gravitational 
plane wave. Exact equations are derived for the observed spectrum, sky location, 
angular size, and brightness of each such source. Some (though not all) of these 
observables have been considered previously in [3], but in forms which were difficult 
to interpret and which could not be immediately compared to known perturbative 
results. Where overlap exists, the new equations obtained here are far simpler both 
to understand and to apply. Throughout, gravitational waveforms and polarization 
content are unconstrained, and light is assumed to follow the laws of geometric optics. 

The basic system we consider is illustrated in Figure 2a. There, a source is 
abstracted to a timelike worldline S and its observer to another timelike worldline 
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O. Images then correspond to past-directed null geodesics from O to S. Given any 
observation event Po & O, exactly one image is assumed to exist in all cases considered 
here||. The remainder of this paper explores the properties of these images. More 
complicated optical problems (such as the interferometer illustrated in Figure 2b) 
can typically be understood as multiple copies of the basic system shown in Figure 
2a. Unless otherwise indicated, the results of this section place no restriction on the 
motions of the source or the observer. 

3.1. Time of flight 

Observations such as those illustrated in Figure 2a provide a natural mapping between 
“observation times” and “emission times.” More precisely, the light rays which connect 
S to O may be used to relate the proper time Tq recorded by O at an observation 
event Po G O to the proper time recorded by S at the emission event Ps G S. 
Such relations play an important role in all optical calculations performed below. The 
redshift or blueshift can, for example, be computed from dr^/dro. Additionally, the 
time difference Stq (or the phase shift) for the interferometer illustrated in Figure 
2 b follows by successively applying the maps Ts(to) applicable to each pair of optical 
elements. 

As a first step towards computing Ts(to) for the simpler system in Figure 2a, 
consider instead the gravitational wave phase Us{uo) at Ps as a function of the 
gravitational wave phase Uo at po. Although the plane wave phase coordinate u is null, 
the difference Uo — Us{uo) may nevertheless be interpreted as a “time of flight” for a 
photon traveling from S to O. This interpretation is supported by noting that i) the u 
coordinate is unique up to a constant affine transformation, and ii) the restriction of u 
to any timelike geodesic is, up to a positive constant, a proper time for that geodesic. 

The phase relation Us{uo) can be computed using Synge’s function a[p,p'). This 
is a two-point scalar which takes two events as arguments and returns one half of the 
squared geodesic distance between those events [25, 29, 30]. Plane wave spacetimes are 
one of the few examples where this function is known explicitly (although generic post- 
Minkowskian approximations are available [31]), with the Brinkmann form appearing 
in, e.g., [18] and the Rosen form in [30]. Given two events p, p' and their associated 
Rosen coordinates {u,v,x^), (u',u',x'*), 

a{p,p') = -{u-u'){v - v') + - x'f{x - x')f (21) 

where ( 7 “^)“^ denotes the matrix inverse of the average ( 7 “^) defined by (15), except 
that this average is to be evaluated between u' and u instead of between uq and u. For 
the trivial case of a spacetime which is globally-flat, there exist coordinates where 'jij = 
and so (Tflat(p,p') = -{u-u'){v-v') + l\x-x'\‘^ = ^[-{t-t')'^+ \x-x'\'^+ {z-z')'^]. 
The general expression (21) for the geodesic distance in an arbitrary plane wave differs 
from this only via the replacement I -G in the transverse directions. 

The existence of the null geodesic which connects the observation and emission 
events illustrated in Figure 2a implies that 

<^{Po,Ps) = criPs,Po) = 0, ( 22 ) 

II There exist special source-observer configurations where pointlike sources appear as extended 
images. Plane wave spacetimes can also admit multiple discrete images on very large scales which 
cannot be described by the single Rosen coordinate patch assumed here [3, 12, 32]. 
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so 

Uo-Us = ^{Vo - - X^fiXo - XsY (23) 

whenever Vs ^ Vo- Using the t and z coordinates defined by (14) instead results in the 
equivalent 

to-ts = ^J[{'y~'^)~^]ijixo - XsYixo - x^y + (zo - Zs)^. (24) 

Supplementing these equations with appropriate parametrizations for the source 
coordinates in terms of the gravitational wave phase allows explicit solutions to be 
obtained for Us{uo) or ts{to)- If the source and observer move on geodesics, the 
necessary parametrizations are given by (16) and (17). 

Once Us(uo) is known, it may be used to relate proper times associated with 
the observation and emission events. The observer’s proper time Tq is related to the 
gravitational wave phase via 

where denotes the observer’s 4-velocity and 1°" = d/dv the gravitational wave’s 
propagation direction. Noting that is Killing, this rate is a constant for geodesic 
observers. More generally, we parametrize it by 

Ao = > 0, (26) 


which coincides with the A appearing in (17) if the motion is geodesic. Regardless of 
acceleration, a (not necessarily constant) Ag can be defined analogously for the source, 
in which case proper times along S and O are related via 


drg 

ATo 




(27) 


If the source and the observer both move on geodesics, Ag/Ao is a constant and this 
equation is trivially integrated to obtain Ts{to) in terms of Us(uo)- 


3.2. Frequency shifts 


The derivative drg/dro appearing in (27) is directly observable. If a physical process 
on S occurs with a characteristic frequency Wg—perhaps the frequency of a spectral 
line or the angular velocity of a pulsar—which is much larger than any frequencies 
associated with the gravitational wave or with the acceleration of S, such a process 
would be observed on O to occur at frequency ujo = (dTg/dro)(jJg. This reduces to 
the familiar Doppler effect in fiat spacetime, but more generally includes curvature 
corrections as well. One way to compute ujo is to implicitly differentiate (23) and 
then substitute the result into (27). Alternatively, note that the gradient of Synge’s 
function <j(jpo,Ps) is tangent to the light ray which connects Ps to Po. It is also parallel- 
transported along that ray, implying that 


Wo _ ^o VgO- 
Wg U^'^bcr 


(28) 


in terms of the observer and source 4-velocities 77“ j. The gradient in the numerator 
here is understood to be evaluated at Po while the gradient in the denominator is 
understood to instead be evaluated at Ps. It follows from (28) that uJo = Wg if 77“ is 
equal to 77“ parallel-transported along the relevant light ray, a generalization of the 
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flat-spacetime result that the Doppler effect vanishes for comoving objects. Such a 
condition is, however, physically unnatural in curved spacetimes; the frequency shift 
is generically nonzero. 

Equation (28) may be directly evaluated using (21) and (26). The result is 
conveniently expressed in terms of the 2-vectors 

V2Xo{uo-u,) ’ V2K{uo-Us) ’ ^ ^ 

where 


Xo^s = Xo - {Uo - Us){j ^)‘y{Uo)Xo, Xs^o = Xs + {Uo-Us){‘y ^)‘y{us)xs. (30) 


All averages here are to be evaluated between Ug and Uo, while, e.g., Xq denotes a 
u-derivative of the observer’s transverse coordinates evaluated at po- Up to overall 
factors included for later convenience, and kg are essentially transverse separations 
between the source and the observer evaluated either at the emission phase Ug or at 
the observation phase Uo- Causality does not allow any observation at u = ito to 
depend on properties of the source at that phase (except for special alignments where 
Ug = Uo), so the relevant observables instead involve an extrapolation of the source’s 
position from u = Ug to u = Uq performed using the geodesic which is tangent to S 
at pg. Comparing (16) and (30), Xg^o represents precisely this kind of “osculating 
extrapolation.” Similarly, Xg^g extrapolates the transverse location of the observer 
from u = Uo to u = Ug using the geodesic tangent to O at Po- In terms of kg and kg, 
the frequency ratio in the presence of an arbitrary plane wave is 


ojg _ dr, _ , 'fi/Mklki - 'f.X{uo)kl,ki\ 

Ws dTo y l + 'J^/{Uo)klki ) 


(31) 


This is valid for arbitrarily-moving sources and observers. It reduces in the flat limit 
jij —>■ 5ij to an expression for the ordinary Doppler shift. More generally, there 
is a sense in which the effects of relative motion and spacetime curvature can be 
disentangled by considering the time derivative of Wo/w, [10, 11, 33]. 

Equation (31) simplifies significantly if the source and observer are assumed to 
move on geodesics which remain at fixed spatial coordinates. When this occurs, 
Ao = A, = 1 and 


kg - kg 


(7 ^( Xg - Xg ) 
\^{Ug Ug ) 


(32) 


The frequency shift then depends on the difference between inverse Rosen waveforms 
at the emission and observation events: 


LOg - UJg 


h^/{ug)-'yij\ug)]klki 


l+7fci {uo)k^gki 


(33) 


One deficiency with this formula is that kg has no immediate physical interpretation. 
It is closely related, however, to the observed position of the source. As explained in 
Section 3.3, that position is naturally described in terms of a unit 3-vector (fc_L,fc||). 
Using (39) to relate this vector to kg, the frequency shift (33) can be rewritten as 

——— = + cos9)~'^k\[I - ^{uo)'y~^{ug)^'^{uo)]k±, (34) 

where 9 denotes the observed angle between the source and the gravitational wave 
propagation direction [cf. (41)]. Additionally, is a square root of in the sense 
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of (13). Although such roots are not unique, a particular choice must be made to in 
order to define an observer rest frame with respect to which the components of k± 
are defined. The appearing in (34) is the same as that used in (38) to construct 
this frame. 

Although exact, (34) is very similar to the first-order perturbative formula 
typically used when discussing gravitational wave measurements via pulsar timing 
[2, 14] or Doppler velocimetry [34, 35, 36]. In terms of the first-order metric 
perturbation ehij which appears in (3), it is well-known that 

——— = :^(1 + cos0)“^[/iy (us) - hij{uo)]k]_yj_ + O(e^). (35) 

UJo I 

The physical meanings of the 0, fcj_, and which appear here are unchanged in the 
exact result (34), although their time-dependence is no longer trivial and relations to 
coordinate quantities are different. The metric difference e[hij(ug) — hij(uo)] appearing 
in the first-order result is, however, generalized to 

- S.^k{Uo)lkl{Us)£,jl{Uo) = [ikHuo) - lZli'^s)%k{Uo)£,0l{Uo)- (36) 
Direct approximations to (34) are discussed further in Section 4.2.3. 

3.3. Source locations 

Gravitational waves may affect not only a source’s apparent spectrum, but also its 
location on the sky. Such locations can be compactly described in terms of a 3-vector 
which resides in the observer’s instantaneous rest frame. More precisely, consider 
an orthonormal triad (e“,e||) which is orthogonal to the observer’s 4-velocity U^. 
It is convenient to align this triad such that one of its components is locked to 
the direction of propagation associated with the background plane wave. Projecting 
the null propagation direction into the observer’s rest frame and normalizing, the 
longitudinal frame vector is then 

e[[ = y2Ao(g“'’ + KU^b = - K- (37) 

It is also convenient to define the remaining two frame vectors by 

el = I"'^(mo) {i{uo)Xo^ + , (38) 

where ^ is a particular matrix satisfying (13). It is easily verified that the resulting 
triad is indeed orthonormal and orthogonal to U^. It is also parallel-transported for 
geodesic observers. If O is not a geodesic, another frame—perhaps one which is Fermi- 
Walker transported—might be more natural. We nevertheless apply these definitions 
in all cases. 

Recalling that the direction of the light ray connecting po to Ps is given by 
—'^a<^iPo,Ps), its rest-frame components describe the observed location of S. If the 
unit-normalized versions of these components are denoted by (fcl,/C||), use of (21) 
results in 

L _ 2^ (uo)ko ~ ^ 1 — fcpq' {uo)ko 

1 + kl'^-^{uo)ko’ " 1 + kl'y-^{uo)ko’ 

where fco is defined by (29). By construction, \k±\^ These expressions hold 

for arbitrarily-moving sources and observers in arbitrary plane wave spacetimes. 
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It is often sufficient to consider only the transverse 2-vector k±. This is 
proportional to ^~^(uo)fco, the first factor of which takes into account that the 
transverse 2-metric 7 (mo) = l^(uo)^(rto) is not Euclidean at the observer’s location. 
The vector kg is more interesting, however. Recalling its definition, the observer’s 
transverse location Xo does not enter on its own, but instead appears via the osculating 
projection Xo^s- Objects therefore appear to be not where they “are,” but where 
they “would have been.” This can be made more clear in terms of the distance 
measurements considered in Section 3.4. Using (42) and the area distance rarea (43) 
in (39), 


fc, = 


det (7 


1/2 


^det7-i(Mo)7-i(us) 


I '^(Wo)(7 ) (Xs-Xo^s)- 


(40) 


It is sometimes convenient to parametrize source locations by angles instead of 
vectors. Let (0, (j)) be polar coordinates constructed in the observer’s sky such that 6 
corresponds to the apparent angle between the source and gravitational wave, while 
(j) = 0 coincides with the direction of (e“ )i. Then, 


fci = 


cos(/)sin 0 \ 
sinsin 0 J ’ 


cosO. 


(41) 


Combining this with (39) shows that the observed latitude 9 explicitly satisfies 


tan^( 6 »/ 2 ) = klj '^{uo)ko. 


(42) 


One interesting characteristic of position measurements is that they can be used to 
deduce the presence of a gravitational wave even when observing very nearby sources. 
While knowledge of k± at any one instant isn’t particularly meaningful, its time 
evolution is. This evolution can be nontrivial no matter how close S happens to 
be to O. It is evident from the presence of 7 “^(uo) — 7 ~^(ms) in (33) that useful 
timing measurements are instead restricted to source-observer separations over which 
significant differences can be expected in the gravitational waveform. 


3.4- Distances 


The next observable we consider is the area distance rarea- If a source’s angular size is 
resolvable and found to be equal to the small solid angle U in the observer’s sky, there 
is a sense in which its physical area is Urarea [12]- This is somewhat technical to derive, 
so we defer to the result obtained in [3] using Brinkmann coordinates. Translating 
that into Rosen coordinates while using (42) results in 

_ V2Xo{uo - Us) I det(7~^) _ 

I -f cos 6 * y a/ det ■y~'^{uo)'y~^{us) 



The first fraction here is essentially an affine distance to the source, while the second 
corrects this by taking into account the expansion of a thin bundle of light rays which 
converge on p^. That latter correction is remarkably simple in terms of the Rosen 
waveform, depending only on the arithmetic and geometric averages of 7 “^ between 
the source and the observer. 

It can also be useful to describe a source in terms of its luminosity distance rium, 
which is related to the area distance via 

drs 
dTo 


rium = 
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One factor of drs/dro occurs here due to the frequency shift experienced by light 
traveling from S to O, while the second arises from considering bundles of light rays 
which converge on ps rather than po. An explicit formula for the luminosity distance 
follows immediately by substituting (31) and (43) into (44). For the special case 
where both sources and observers move on geodesics at fixed spatial coordinates, the 
frequency shift is given by (34) and hence 



(45) 


4. Perturbative approach 

Equations (23), (31), (40), (43), and (44) provide exact prescriptions for the time 
delays, frequency shifts, positions, and distances of generic sources in plane wave 
spacetimes. We now consider their perturbative expansions. This is done for two 
reasons: First, it provides a clear connection between the exact results derived here 
and the various approximations which have appeared in literature. Second, some 
physical implications of the optical formulae are more easily understood when written 
in an approximate form. In particular, the results of Section 3 can be expanded beyond 
the lowest-order approximation which has typically been considered in the past, and 
doing so demonstrates that some higher-order corrections grow very rapidly with the 
source-observer distance. 

Two types of nonlinearity arise here. The first of these has its origin in the 
nonlinearity of Einstein’s equation as applied to the Rosen-type plane wave metric 
(12) and its associated waveform "fij. Perturbative expansions of this waveform are 
considered in Section 4.1, and are sufficient to qualitatively understand that nonlinear 
effects can accumulate at large distances. The various gravitational wave observables 
are, however, nonlinear functionals of 'jij. This is taken into account when deriving 
explicit second-order expansions for those observables in Section 4.2. 

4-1. Expanding the metric 

As explained in Section 2.2, gravitational plane waves in general relativity can be 
trivially described in terms of the Brinkmann waveform Hij which appears in the line 
element (6). This directly determines the curvature, and is restricted by the vacuum 
Einstein equation only to be trace-free. The observables discussed in Section 3 are 
not, however, written in terms of Hij. They instead depend on the Rosen waveform 
7 y , which is related to Hij via the nonlinear and nonlocal expressions (11) and (13). 
Although the vacuum Einstein equation could be applied to find a nonlinear differential 
equation for the Rosen waveform alone, we take the perspective that it is more natural 
to instead describe a plane wave in terms of Hij and then to derive jij from that. 
It is also convenient to view the perturbative expansion performed here not as an 
approximation for a single system which involves a small quantity e, but rather as an 
approximation for an entire family of systems which are smoothly parametrized by 
e. Orders in perturbation theory then correspond to differentiations with respect to e 
evaluated at e = 0. 
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The precise family of plane wave spacetimes considered here^ is defined to be the 
set of Brinkmann line elements (6) with 

(u; e) = eSjij(u), e > 0. (46) 

The “reference waveform” Sjij{u) = d^Hij(u] e) is assumed to be fixed and independent 
of e. Additionally supposing that tri 3 (u) = 0, it follows from (7) that each spacetime 
in this family is an exact solution to the vacuum Einstein equation. The parameter e 
controls the amplitude of a wave’s curvature, but not its polarization, frequency, or any 
other local properties. As expected from such an interpretation, different members 
of this family are physically distinct: Choosing any nonzero e and any e' ^ e, the 
waveform Hij(u]e') = {e'/e)Hij{u]e) cannot be transformed into Hij(u\e) using the 
gauge transformation (9). It is also clear that the wave disappears entirely entirely in 
the e —>■ 0 limit. 

Given a trace-free Sjij{u), an associated Rosen waveform jij(u;e) may be found 
by first constructing a matrix e) with the expansion 

4(m; e) = 4(0) (w) -I- ... -I- e”4(n)(if) +- (47) 

Applying (11), 

4(0) = 0, 4(n) = ■Q4(n-1) (48) 

for all n > 1. Iteratively solving these equations then results in an approximation for 
^ij{u;e). Using it together with (13) produces a family of Rosen waveforms 
with the perturbative expansion coefficients 

7(n)(u) = ^ i9,"7(u;e)|^^o = ^ 4 Jp){u)4(9)(m). (49) 

p+q=n 

Recalling (3), we use the more standard notation h{u) interchangeably with 7(i)(m) 
for the first-order metric perturbation. 

As already emphasized in Section 2.3, different Rosen waveforms—corresponding 
to different initial conditions for (11) —can be used to describe the same physics. 
Ambiguities of this kind are easily resolved for the sandwich waves described in 
Section 2.4, in which case it is natural to fix a particular waveform by demanding 
that = 4ii = Sij in one of the locally-flat regions. More generally, however, there 
does not appear to be any “preferred” choice. 

Our approach is to fix a convenient solution for ^ij, and then to note that all 
other possibilities are related via (19). Letting this particular solution satisfy the 
non-perturbative initial condition 

^ij{uo;e) = Sij, ^ij{uo;e)=0 (50) 

for some constant wq, it follows from (48) that ^(o) = I and 

pU pUi 

4(n)(w) = / dui du2Sj{u2)^(n-l){u2) (51) 

Juq J Uq 

and for all n > 1. Hence, 

tr4(n) = 4^„) - 4(n) =0 (n odd), (52) 

% It can sometimes be interesting to instead consider families of waves where the curvature depends 
smoothly on e, u and u/e. Any portions involving u/e vary rapidly as e ^ 0, thus evoking the 
concept of a “high-frequency limit” [37, 38]. Considering such families is equivalent to a type a 
singular perturbation theory [as opposed to the regular perturbation associated with (46)]. See also 
the last paragraph of Section 4.1.2. 
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and 


^(") + ^(n) = (tr€(«))-? (n even). (53) 

Combining these constraints with (49) shows that each 7 („) generated by a ^ which 
satisfies (50) is symmetric and that 

tr7(n)=0 (nodd), 7 („) = ^(tr 7 („))J (n even). (54) 

The familiar tracelessness of the first-order TT-gauge perturbation h = 77 ) therefore 
generalizes to all odd-order metric perturbations. It does not generalize to even orders; 
those expansion coefficients are instead “pure trace.” 

Through second order, these results imply that the particular Rosen waveform 
defined by (50) is explicitly 


7 (m; e) = I -I- 2e f dui 
Juo 


du 2 Sj{u 2 ) + tr 


•1 pU\ 

dwi / Au2^{u2) 

1 j Uc\ 


2 


pU pUi pU 2 pU 3 

+ / dui / dM2 / dW3 / dU4Sj{u2)Sj{u4) 


/ + 0 (e^). 


(55) 


The second line of this equation generically grows with u, which can be seen more 
clearly if the waveform is rewritten in terms of the first-order perturbation /17 instead 
of S)ij = ^hij. Making this replacement and integrating by parts results in 


7 (m; e) = 7 -I- eh{u) + -e^ tr 


h^{u) — / dwi / du2h^{u2) 

J Uq ^ Uq 


I+ 0(6^). 


(56) 


The trace of cannot be negative, so the double integral here—and therefore the 
second-order metric perturbation as a whole—typically grows with u. Indeed, there 
exists an e-dependent lengthscale beyond which the second-order term outpaces the 
first. We now consider two examples which demonstrate this explicitly. 


4-1.1. Gravitational wave bursts with memory Consider a sandwich wave as described 
in Section 2.4. More specifically, consider an e-dependent family of waves (46) with 
^ij{u) zero everywhere except in a small region u S (—(5,0) inside of which the 
curvature oscillates several times. Choosing uq < —6 in (50) guarantees that the 
metric before the burst is trivial: ^ij(u\e) = 6ij. The metric after the burst must, 
however, be generated by the square of ( 20 ), where aj^^{e) and bjr^{€) are constant 
matrices determined by the details of the curved region. If ajr^ (e) ^ I or bjr^ (e) ^ 0, 
the gravitational wave burst is said to exhibit memory [39, 40] in the sense that it 
permanently displaces initially-comoving pairs of freely-falling test particles. It follows 
from (51) that the first and second-order perturbations are 

^ pO pui 

= ^b,{0) = J dui J du2Sj{u2), (57) 

ibr+){i) = ^^( 0 ) = J duiSj{ui), (58) 


{aj^+)(2) 

(^•^+)( 2 ) 


I 


^ pO pui 

(aj^+)(l) ~ 2 J /“l J ^dM2/l^(M2) 

0 


I 


tr 


(aj^+)(i)(bj^+)(i) - T / dwi/r (mi) 


1-6 


I, 

I. 


(59) 

(60) 


and 
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Although gravitational wave bursts with memory can arise from many physical 
scenarios, classic discussions [40, 41] consider those waves which are emitted either 
from the scattering of multiple masses or the explosion of one mass into multiple 
unbound components. Applying the quadrupole formula in such cases suggests that 
(bjr^)(i) = 0.+ Assuming this, a coordinate rotation may always be performed 
together with a rescaling of e such that 

(a.F+)(i) = ^ 0 -1 ) ■ 

Now consider second order effects. It follows from (60) that even though bjr^ vanishes 
at first order, it cannot vanish at second order. Indeed, tr(b_ 7 r^)( 2 ) < 0 in these cases. 
Applying (13) and (20), the second-order waveform is explicitly 

7 (u;e) = /-he|(^p [l + 4 tr(a^^)( 2 )] jj 

+ e\[tr(b^J( 2 )]/ + 0(e3) (62) 

when Uo > 0. If the Hrst-order waveform inside a burst is schematically of the form 
h{u) ^ cosuju and oscillates N ^ ujS times, it follows from (60) that tr(bjr^.)( 2 ) 
—Nuj. The second-order metric perturbation therefore becomes comparable to the 
Hrst-order perturbation when 

oju ^ {eN)~^ 1. (63) 

Of course, 7 (i) and 7 ( 2 ) are qualitatively different when this occurs. One is trace- 
free while the other is a pure trace. On even larger scales where oju ~ (e^A^)“^, the 
determinant of the Rosen waveform goes to zero. Recalling the line element (12), the 
metric itself becomes singular on this scale. Such effects are not an artifact of the 
perturbative expansion, but instead signal the breakdown of the Rosen coordinate 
system due to the gravitational focusing of nearby geodesics. 

4-1.2. Monochromatic waves Another instructive example is provided by a 
linearly-polarized gravitational wave with constant amplitude. Let this wave be 
monochromatic with angular frequency oj in the sense that it can be described by 
the family (46) with 

■5(w) = Y ( Q ^ cos OJU. (64) 

It follows from (11) and (46) that if a; = 1, the 2, 2-component of ^ij(u;e) satisfies 
the Mathieu equation ( 1 ). [The equation for general oj differs from this only via 
the rescaling u —>■ wu, resulting in ^ 22 {u\e) + \eoj^ ^ 22 {u] e) = 0.] We now show 
that important features of the metric perturbations associated with monochromatic 
gravitational waves are captured by the approximate Mathieu solution (2). 

Although it is possible to express all components of jij (ir, e) exactly in terms 
of Mathieu functions, consider instead the perturbative expansion described above. 
This first requires imposing initial conditions for ^ij{u',e). Unfortunately, there is no 
natural phase at which to apply (50). Indeed, all choices are more complicated than 

+ The condition = 0 implies that through first order, pairs of test particles which are 

comoving before the burst are also comoving after the burst; only their displacements might be 
permanently affected at 0(e). Some cases where an additional “velocity memory” arises are discussed 
in [42]. See also [41, 43] for an electromagnetic memory effect where nontrivial velocity changes occur 
even at lowest order in simple scattering problems. 
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Figure 2: Plot of 722(w) for a monochromatic gravitational wave with S^ij(u) given 
by (64) and e = 10“^. An exact solution is shown together with its first and second- 
order approximations. Initial conditions for all three curves are matched at u = 0. 
The 0{e) approximation is poor after = 10 curvature oscillations, although the 
O(e^) approximation still works well there. A coordinate singularity appears in fewer 
than e~^ = 100 oscillations. 


necessary. Setting uq = 0 for definiteness, (49) and (51) imply that hij involves an 
oscillating term as well as a constant offset. Eliminating this offset by modifying the 
initial condition at 0(e) recovers the expected first-order perturbation 

h{u) = ~ ^ Q ^ cos uju. (65) 

The second-order perturbation is then 

7(2)(u) = ^[l-2(wu)^-b3cos2a;M]/-b2[|(2)(0)-fui(2)(0)], (66) 

where the final two terms allow for O(e^) modifications to the initial condition (50). It 
is clear, however, that no matter what these modifications are, the quadratic growth 
when |a;u| ^ 1 cannot be eliminated. This is a second-order effect which becomes 
comparable to first order effects when 

\lou\ ~ £-1/2 > 1. (67) 

It can be understood intuitively as a version of (63) where the number of gravitational 
wave oscillations N isn’t fixed, but is instead of order uju. For a given wave amplitude 
e, it follows that nonlinear effects tend to be more important for approximately- 
monochromatic waves than for bursts. Also note that as suggested by (66), Rosen 
coordinates break down when |a;rt| ^ e~^ >> £“1/2. This is illustrated explicitly in 
Figure 2 for the special case e = IO-2. 

The large nonlinear effects considered thus far do not oscillate. Monochromatic 
waves where £ <C 1 first experience relatively-large oscillatory corrections at 0(£^). 
Assuming initial conditions for ^(i) such that (65) holds, the large-u portions of 7(2) 
and 7(3) can be viewed as conformal corrections to the first-order waveform: Using 
an ellipsis to denote omitted terms at second and third order which grow more slowly 
than u^, 

^ cos UJU 


+ ... + 0{e^). (68) 


7(m;£) = 


1 — -{eujuY 


I- 


1 0 
0 - 
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The growing nonlinearities discussed here can be interpreted in some ways as 
manifestations of an averaging effect which occurs more explicitly in “high-frequency” 
perturbations where families of waveforms e) are considered which are smooth 

in u/e and e, but not in u and e. In that context, the zeroth-order metric is no 
longer Minkowski, but instead satisfies the non-vacuum Einstein equation with an 
effective stress-energy tensor sourced by the gravitational wave [37, 38]. Indeed, 
this stress-energy tensor is controlled by an average of the same “energy density” 
which is responsible for the second-order metric (56) in an ordinary perturbative 
expansion. The high-frequency viewpoint is not pursued here firstly because its 
formulation requires choosing one small parameter which simultaneously controls both 
wave amplitudes and distances. This somewhat limits flexibility and clarity for the 
questions considered here. More importantly, different observables can react very 
differently when considering non-smooth families of spacetimes. The usual high- 
frequency formulations in general relativity are designed to compute metrics which 
are limits of solutions to Einstein’s equation. It is not necessarily true, however, that 
limits of optical observables can be easily described in terms of these same metrics (or 
even that their limits exist at all). 

4-2. Gravitational wave observables: perturbative expressions 

Perturbative expansions for the various observables computed in Section 3 may now be 
found for families of plane waves described by (46). Completely general expressions 
would be quite complicated, so we restrict attention only to those cases where the 
observer O remains fixed at the spatial coordinates {xo,Zo)- Also suppose that the 
source S remains fixed at {Xs,Zs). As guaranteed by the construction of the Rosen 
coordinate system described in Section 2.3, all such worldlines are timelike geodesics. 
They are not preserved, however, by changes in the Rosen waveform; coordinate 
transformations (19) which modify 7^ also impart initially-stationary geodesics with 
nonzero coordinate velocities. Our assumption that particular sources and observers 
remain at fixed spatial coordinates is therefore coupled to the choice of a particular 
Rosen waveform. Although no detailed prescription is used here, we do assume that 



These equations are implied by the initial condition (50), although they also hold more 
generally. 

4-2.1. Zeroth order Our first step is to compute all observables in the flat-spacetime 
limit e —> 0. Denoting orders in perturbation theory by analogy with (47), the zeroth- 
order observables are 



(70) 

(71) 

(72) 

(73) 


As expected, both the area and luminosity distances reduce to the coordinate distance 
r at this order. Frequency shifts also vanish, and a source’s apparent location in the 
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observer’s sky is given by the expected function of its coordinate position. Somewhat 
less obviously, (23) implies that the zeroth-order emission phase is 

= Wo - ;^(1+ cos0(o)) (74) 

in terms of the observation phase Uo- 

All higher-order perturbative expressions obtained below depend in various 
ways on these zeroth-order quantities. Consistently retaining “(0)” subscripts or 
superscripts on each of them would considerably increase clutter, so this notation 
is suppressed when no confusion should arise. 


4-2.2. Time of flight The first-order perturbation = 9eUs|e=o to the zeroth- 
order emission phase (74) is easily found by differentiating (23) with respect to e while 
holding Uo fixed. This results in 


dr = 


(75) 


where k± is understood to refer to the zeroth-order source direction Similarly, 

(hij) denotes an average of the first-order metric perturbation between and 
Uo- Decomposing hij into its -I- and x components using (4), the emission phase 
perturbation can alternatively be written as 




r 


1 

2a/2 


[(/i_|_) cos 2(1) + (hx) sin 2(j)] sin^ 9. 


(76) 


This vanishes for all sources which are aligned (0(o) = 0) or anti-aligned (0(o) = tt) 
with the gravitational wave. If that wave is linearly-polarized, coordinates may be 
chosen such that hx = 0, implying that also vanishes for any source which is 
nominally located on one of the four meridians 


<(’(0) = (1 + 2n)7r/4, n = 0,1,2,3. 


(77) 


Regardless of polarization, the averages appearing in (76) are typically small if the 
waveform is approximately oscillatory and there are many oscillations between the 
source and the observer. These averages can be important, however, when there is 
significant gravitational memory. 

The second-order perturbation = ^d^Usle=o to the emission phase is 

somewhat more complicated. Differentiating (23) a second time with respect to e 
shows that 


( 2 ) 

Us 


r 


^tr [(7(2)) + (/i)"- 

V-s f hijiUs) {^ij) 
2r \ 1-1- cos 9 


{h'^)] sin^ 9 


(78) 


( 2 ) 

Unlike its first-order analog, it) ’ does not tend to zero for oscillatory, memory-free 
waves at large distances. Part of it is also independent of . 
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4 . 2 . 3 . Frequency shifts A perturbative expansion for the relative frequency shift 
(wo — u}s)/oJs induced by a gravitational plane wave may be found by differentiating 
the exact result (34) with respect to e while holding Wg fixed. At first order, 


^( 1 ) 1 

= ;:(1 + COS0) ^[hij{us) - htj{uo)]k]_k]_, 
UJs I 


(79) 


which is immediately seen to be equivalent to the well-known [2, 14] approximation 
(35) when combined with (71). In terms of /i+ and /ix) h is rnore explicitly 

^(1) 1 

-2— = -(1 - cos 9) {[h+{us) - h+{uo)] cos 2 (;ii -|- [hxius) - hy^{uo)] sin 2 (()}. (80) 

CJe 2 


Equation (74) implies that — )■ Uo as d(o) — t tt, so uj'lf' = 0 for any sources which 

are aligned or anti-aligned with the background gravitational wave. If the wave is 
linearly-polarized and coordinates are chosen such that hx = 0 , first-order frequency 
shifts also vanish when ^(o) satisfies (77). 

Regardless of polarization, differentiating (34) again with respect to e results in 
the second-order correction 

^( 2 ) 2 

= t(1 - cos 6 i)tr { 7 ( 2 )(ms) - 7 ( 2 )(uo) + [h{us) - h{uo)][2{h) - h{us) - h{uo)]} 
LOg 4 


,( 1 ) 


^(1-|-COS0) ^ {uj^^'>/ujs)[hij{us) - {3 + cos9){h^j)]+ui'^'>h^j{us) 


k±k^±. (81) 


The second-order metric difference 7 ( 2 ) (^1°^) — 7 ( 2 )(mo) which appears here is directly 
analogous to the first-order difference — h{uo) found in Wo^^/wg. Also note that 

the final term involving h(us) may be interpreted as a straightforward correction to 
the first-order expression obtained by the replacement -I- ...). 

Overall, and both vanish when 0(o) —>■ 0 or 0(o) —)• tt. 


4 . 2 . 4 . Source locations The apparent location of a source is governed via (41) by the 
unit vector (fcj_,A:||) introduced in Section 3.3. Differentiating the exact expression 
(39) for with respect to e shows that first-order astrometric effects associated with 
a gravitational plane wave follow from 


= [(h) - ^Huo)]k± + - 


hij{Uo) i^ij) 
1 + COS 9 


(hij) 




fc_L, 


(82) 


a result which has also been derived (using different methods) in [4]. One of its 
consequences is that the angle 9 between the source and the gravitational wave is 
perturbed by 

9(1) = [(hy) - hij{uo)] + {hij) cosdjfc^L^i CSC 9. (83) 

If a wave is linearly-polarized with hx = 0, sources which nominally lie on the 
meridians (77) experience vanishing latitudinal motion. These sources can, however, 
appear to rotate slightly around the gravitational wave propagation direction: 9(i) = 0 
but (j)(i) = ±[{h+) - i/i+(Mo)]. 

Regardless of polarization, two special cases of (82) may be understood 
immediately. The first of these supposes that (hij) can be neglected, and may be 
thought of as a “large-r” limit in the presence of memory-free waves. Applying it, 

—)■ ~ COS0) [hj.{uo) cos2(/) + hx (uo) sin2(()] J — h{uo)'^k±. (84) 







Optics in a nonlinear gravitational plane wave 


24 


Gravitational waves can also produce nontrivial astrometric effects in a “small-r” limit 
where (hij) —>■ hij{uo) and 

> i|h,(uo) — sin^ 9\h+{uo) cos 2(j) + hx (uo) sin2(/)] (85) 

In either of these cases, first-order position perturbations depend on the waveform only 
at the observer (where u = Uo). The angular dependence of this effect is nevertheless 
distinct for near and distant sources. 

If a gravitational wave is linearly-polarized, the matrices which multiply in 
(84) and (85) depend on Uo only via an overall scaling. Sources in either limit therefore 
appear to move coherently along straight lines in the presence of linearly-polarized 
gravitational waves. The orientations and relative amplitudes of these lines depend, 
however, on each source’s nominal location (0(o),0(o)). More complicated apparent 
motions can arise in either the small-r or large-r limits for gravitational waves which 
are not linearly-polarized, and even for linearly-polarized waves when {hij) is neither 
negligible nor approximately equal to hij{uo). 

Astrometric effects are significantly more complicated at higher orders. For 
brevity, we therefore present second-order corrections only for the latitude 9. 
Expanding (42), 


5'(2) = ^tr|[(7(2))-7(2 )(uo)] -f (7(2))cos6»-h (I+ cos6 I) [{h)^ - {h^)] 


[(/i) — h(uo)]^ sin0 


- hij(us) 


V2'. 


k±k]_ — -{2 — cos 9 )9^1^ > csc 9. 




( 86 ) 


All terms involving 7(2) in this equation are easily seen to be direct generalizations of 
the first-order metric perturbations appearing in ^7) [cf. (83)]. 


4-2.5. Distances The remaining observables considered here are the area distance 
Tarea and the luminosity distance num- First-order perturbations for these quantities 
may be found by expanding the exact expressions (43) and (44), which results in 


^(1) 

' area 


— 2 ( ^^^3) + 


{hij) hij{uo) 

1 + cos 9 


P 


(87) 


and 

( 1 ) ( 1 ) 

+ {1 +cos9)~^[hij{uo) - hij{us)]k]_yj_. ( 88 ) 

The second of these results describes how an object’s apparent brightness is affected 
by a gravitational wave. It has sometimes been stated in the literature that gravity 
first affects brightnesses only at second order in general relativity (e.g., [5, 44]). A 
typical argument appeals to the Raychaudhuri equation, which can be used to show 
that the expansion of a null congruence is unperturbed through first order for objects 
at a fixed affine distance. This is misleading, however. Gravitational waves (and more 
general geometries) do affect affine distances at 0(e). This and the first-order time 
dilation both contribute nontrivial first-order perturbations to Carea and rium- 

Second-order distance perturbations must take into account changes in a source’s 
affine distance, time dilation, and, unlike in the first-order case, the gravitational 
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focusing of null congruences. All of these effects are taken into account automatically 
by differentiating (43) twice with respect to e: 


„(2) 

' area 


= ;^tr|7(2)(uo)cos6» + 7(2)(Ms) - (1 + cos 0) cos 0 [(7(2)) + {h)^ - {h'^)] 
+ (1 - cos9)[{h) - hiuo)]"^ - ^[h^{uo) + - 2{h)‘^]j + 


1 


1 + COS 6 


+ -(3 + cos9){hij) - hij{us) 


P 


(89) 


The second-order perturbation to the luminosity distance follows from this, (79), (81), 
and (87) via 


,( 2 ) 

lum 


( 2 ) ( 1 ) 

^ - 2(a;(^)M). (90) 


4-3. Summary of second-order effects 


Most terms in the second-order expressions derived above are relatively uninteresting 
in the sense that their magnitudes are comparable to squares of typical first-order 
magnitudes. If these latter 0(e) effects are only marginally detectible, their squares 
are hopelessly small. More interesting are the second-order terms which can acquire 
large numerical coefficients. As motivated in Section 4.1, 7(2) can be large, thus 
implying that (78) and (86) simplify to 


u?) = - 


^tr(7(2))sin2 


-H .... 


8 ( 2 ) = 4 tr{[(7(2))-7(2 )(uo)] + ( 7 ( 2 )) cos 0} sin 0-1-..., 


(91) 

(92) 


at large distances. These are identical to the expressions which would be obtained if the 
first-order expressions (75) and (83) for and 0(i) were applied with the substitution 
h —> £7(2); providing a sense in which the potentially-significant contributions to Us 
and 0 through second order mimic first-order effects with an “effective first-order 
metric” h £7(2). This effective metric has a nonzero trace, and therefore mimics a 
third type of gravitational wave polarization—sometimes referred to as a “breathing 
mode” [45] —which introduces a distinct ())-dependence into various observables. It 
must be emphasized, however, that this “extra polarization” is only an effective 
phenomenon. It is entirely determined by the ordinary -|- and x polarization states, 
and therefore does not represent a physically-independent degree of freedom. This 
is also an effect which arises only when considering nonlocal observables, and not in, 
e.g., direct local measurements of Rabc^- 

Although the notion of an effective metric applied to first-order expressions 
describes all potentially-large second-order contributions to and 0, it cannot do 
so for all interesting terms 


,,(2) 1 1 

= -(l-cos0)tr[7(2)(u^)-7(2)(Mo)] + -(u(iVO 


rhtj{us)k]_k]_ 


LO 


1 -I- cos i 


+ ... 


(93) 


associated with the second-order frequency perturbation (81). The trace terms here 
are indeed those which would be obtained by adding an appropriate correction to 
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the metric perturbation appearing in The remaining portion of the second- 

order frequency perturbation is different, however. It arises from the wave-induced 
perturbation to the emission phase, and can be important when hij varies significantly 
over scales of order Somewhat more precisely, use of (75) shows that the hij 

term in (93) can be large if {hij) is significant and there is a sense in which rhij ^ hij. 
This occurs if, e.g., a gravitational wave simultaneously possesses both high and low 
frequency components. See Section 5.3. 

The effective metric concept fails completely to describe the second-order 
perturbations to the area and luminosity distances. Using (89) and (90), these 
perturbations are dominated by 
( 2 ) 1 

= - tr [7(2 )(uo)cos 6' -b7(2)(ws) - (1 + cos d) cos 0(7(2))] -b ... (94) 

and 

^(2) 1 

= 4 “ cos 0)7(2 ) (mo) + (2 cos 0 - 1)7(2)(Ws) 

- (1-b cos 0) cos 0(7(2))] - 

at large distances, which could not have been guessed from the first-order expressions 
(87) and (88). This is because the gravitational focusing of neighboring null geodesics 
is essential to both of these expressions, but has no first-order analog in vacuum general 
relativity. 

5. Gravitational wave examples 

We now consider the optical effects of three types of gravitational waves in order 
to illustrate some physical consequences of the perturbative expressions derived in 
Section 4. The first of these examples involves a gravitational wave which is both 
monochromatic and linearly-polarized. Next under discussion is a fast burst with 
nontrivial memory. Lastly, we consider a superposition of these two possibilities. 

5.1. Optics in a monochromatic wave 

Perhaps the simplest physically-interesting gravitational wave is a linearly-polarized 
example whose curvature is monochromatic in the sense described in Section 4.1.2. 
Specifically, consider a family (46) of waves with curvatures given by (64). Also 
suppose that Rosen coordinates have been chosen such that 7(1) = ft, is given by (65) 
and 7(0) = /. 

The first-order frequency shift for freely-falling sources and observers embedded 
in such a wave and remaining at fixed spatial coordinates follows from (80): 

= ^{cosuJUo — cosuJUs){l — cos9) cos2(j). (96) 

This generically oscillates as Uo is varied, but vanishes for any sources with ())(o) 
satisfying (77). Other observables depend on the average waveform (hij) between 
the emission and observation events. Using (74) to define the zeroth-order estimate 

^(1 + COS0) (97) 


N = 
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for the number of gravitational wave cycles between these events, (hij) ^ N ^ over 
large distances where ^ 1. It follows that averages can be ignored at first order 
in the large-distance limit. Equation (75) then implies that ui^'^ —> 0. The first-order 
position change is nontrivial, however, and may be described by 

0(1) —>■ - cos wuo sin 0 cos 2((), (/ii(i) —>■ — - cos wMq sin 2(/). (98) 

It also follows from (87) and (88) that the first-order perturbations to a source’s 
apparent distance are 


^arL/^ i(l — COS0) cos2^cosa;uo, (99) 

’’lum/^ —>■ (1 — COS0) cos2(/)(cosa;'Us — ^ cosujUq), (100) 


when (hij) can be ignored. 

Continuing these calculations through second order, it is implied by (66) that the 
dominant contribution to the waveform is 

7(2){u) = -^(a;u)2/-t... (101) 


at large u. Substituting this into (93) results in 


u„ - u. 


(1 — COS0) -I-_ 


( 102 ) 


The factor appearing here can be enormous at large distances, potentially 

allowing the magnitude of to compete with The temporal and angular 

dependencies of the first and second-order effects are very different, however. Similar 
comments also apply to the other observables considered here. 

Consider, for example, 0(2). This involves the second-order average (7(2))) which 
is not generically negligible at any distance. Computing it using (101) and substituting 
the result into (92) shows that 

TT^iV^ 


72) 


12 


- 1 


(1 -I- COS0) 


- 3 


Uo - Us 


COS 0 > sin 0-1- 


(103) 


which is again proportional to Although the first-order observables oscillate 

for monochromatic waves, their second-order counterparts act (at least over short 
observation times) more like offsets. These offsets depend in a characteristic way on 
both the distance to a source and its angular separation from the wave propagation 
direction. Their magnitudes can be comparable to first-order effects when N ^ cor 

e-1/2. 


5.2. Optics and the memory effect 

As a second example, consider a short burst of gravitational waves as described in 
Section 4.1.1. We do not model the burst itself, but only its memory in the form 
of the second-order Rosen waveform (62). Observations may then be split into three 
main phases. These are the i) early times where Ug < —S, ii) intermediate times where 
Uo > 0 but Us < —6, and iii) late times where Uo > 0 and Us > 0. The first and last 
of these phases involve light propagating entirely through flat regions of spacetime. 
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We start by considering first-order effects. The frequency shift is particularly 
simple, vanishing at both early and late times, while holding the constant value 

/ojs =--{I — cos9) cos2(j) (104) 

at intermediate times. The true frequency shift would not, of course, jump between 
these possibilities instantaneously. Transitions would instead last for observation times 
of order the burst time S, and would depend on detailed properties of the waveform. 
Also note that (104) suggests that does not vanish as 0(o) This is in 

conflict with the general comments following (80) , and is an unphysical artifact of the 
discontinuity introduced in the waveform by ignoring timescales of order 6. 

Optical observables other than the frequency shift depend on {hij), which 
cannot necessarily be ignored when memory effects are significant. Indeed, this 
average imparts various observables with the “continuous component” of their time 
dependence. For example, (75) implies that u's'^ is initially zero, changes linearly with 
Uo via 

=- ^ f — — —^ sin^0cos2^ (105) 

2v2 \uo - Us J 

at intermediate times, and then saturates to a constant value at late times. 

Astrometric effects are somewhat more complicated. Through first order, a 
source which is initially stationary at (9,(j)) = (0(o), <)>(o)) rapidly moves to a position 
determined by 

0 ( 1 ) = — i sin 0 cos 2 (/), (/i(i) = i sin 2 (/) (106) 

when Uo = 0. These perturbations then vary linearly with Uo until saturating at 

0 ( 1 ) = i sin20cos2(/), (/)(i) = — ^ sin2^, (107) 

where they remain at late times. That the asymptotic angular perturbations are 
nonzero is a consequence of the finite displacement memory associated with the 
gravitational wave. 

The distance observables rarea and rium are both equal to r at early times. They 
then suffer “immediate” equal and opposite first-order perturbations 

4nL = -^irL = \r{l - COS 0) cos 2(j) (108) 

when Uo = 0. These perturbations subsequently vary linearly with Uo until reaching 

= ^^sin^^cos2</), r|,^^^ = rWa + ?'(l-cos0)cos20 (109) 

when Uo approaches r{l + cos 0(o))/‘\/2 from below. The perturbation to the area 
distance retains this value for all later times, while the perturbation to the luminosity 
distance instead jumps so as to agree (again) with riHa at late times. 

Continuing these calculations through second order, we focus on the contribution 
due to the constant tr( 6 jr^)( 2 ) appearing in (62). As explained following (60), 
this constant must be negative if the first-order memory is to perturb only the 
asymptotic positions of test particles, and not their velocities. The relevant portion 
7 ( 2 ) (w) = it[tr(b_ 7 r^)( 2 )]i' -I- ... of the second-order metric then contributes 

= -^Uo[tr(6^_^)(2)](l - COS0) -f ... 


(110) 
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to the frequency shift at intermediate times. This is nowhere negative, and therefore 
corresponds to a linearly-increasing blueshift. It is an effect which saturates at 

= -^[tr(b^^)(2)]sin2 6 »-b..., (Ill) 

where it remains for all late times. Gravitational wave bursts therefore induce 
permanent blueshifts at this order: Initially-comoving test particles are focused by 
the wave, experiencing a small “velocity memory” at late times. This is a nonlinear 
effect proportional to the nominal source-observer distance r. 

Similarly considering the dominant large-distance effect on 9 at second order, it 
follows from (92) that 0(2) continuously changes from zero according to 


5'(2) = -Uo[tr(6^_^)(2)] 


(1 + cos 9) — 2 


sin0 -I- ... (112) 


at intermediate times. This is nowhere negative, indicating that second-order effects 
tend to make objects appear to bunch up against the wave propagation direction. This 
behavior transitions to 


0(2) = -(uo-u,)[tr(bjrJ(2)l 


cos 9 — (1 -I- COS0) 


sin0-|-... (113) 


at late times, the sign of which depends on 0(o) and Uo- The nontrivial late-time 
dependence of 9 on Uo is another consequence of the velocity memory imparted by the 
gravitational wave burst at second order. 


5.3. Superimposed bursts and continuous waves 

Our last example consists of a superposition of the monochromatic and burst-type 
waves discussed above. Allowing the relative amplitudes of these components to differ 
by a constant factor h^o, suppose that except in a small neighborhood of u = 0, 

h{u) = [/ioo0(u) - coswrt] ^ Q -1 ) ’ 

where 0(u) is the unit step function. All first-order perturbations to the 
optical observables in this case have the form (results of Section 5.1) -|- hoc x 
(results of Section 5.2). 

Second-order effects are potentially more interesting. In particular, the hij 
contributions to the frequency shift (93) and the luminosity distance (95) can be 
significant. Using (105) while assuming that the number of cycles (97) associated 
with the continuous component of the wave satisfies A^ ^ 1, this reduces to 

^ =—7rA^/ioo(l — cos 0)^ cos^ 2()) sin wits (115) 

at late times. Unlike the other second-order examples considered here, this oscillates 
with the same frequency as the first-order waveform. The prefactor irNhao can also 
be extremely large, particularly if the memory amplitude is much larger than the 
amplitude of the monochromatic component so hao S> 1. 


/ rhij{us)k]_yj_ 
r I 1-1- cos 0 
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We have derived the exact time delays (24), frequency shifts (31), sky positions (39), 
area distances (43) , and luminosity distances (44) associated with optical observations 
in the presence of arbitrarily-varying gravitational plane waves. Together, these results 
provide a simple and non-perturbative framework with which to explore the physics 
of nonlinear gravitational waves. 

One conclusion is that the optical effects associated with gravitational plane 
waves appear particularly simple when those waves are described in terms of (u ), a 
nonlocal variable which generalizes the familiar waveforms of TT-gauge perturbation 
theory. Non-perturbatively, is a square root of the transverse metric jij = 
in a Rosen-type coordinate system. It also represents the nontrivial components of a 
Jacobi propagator, and therefore describes separations between neighboring families 
of geodesics. The physical character of ^ij can be understood by noting that it satisfies 
an ordinary differential equation = Hik^kj which also describes a set of coupled 
parametric oscillators. The instantaneous natural frequencies of these oscillators 
directly correspond to those curvature components Hij{u) = which 

represent the freely-specifiable degrees of freedom associated with the gravitational 
wave. Although the equation satisfied by is linear, its solutions depend nonlinearly 
on the “local waveform” . It is through this nonlinearity—which is more closely 
connected to the geodesic equation than to Einstein’s equation—that interesting 
optical effects can arise on large scales. 

Much of our discussion examines these effects perturbatively. An expansion 
+ •. ■ for the metric square root is obtained in Section 
4.1, where e is a small parameter which controls only the overall scale of Hij. The 
various optical effects considered non-perturbatively in Section 3 are then specialized 
in Section 4.2 for freely-falling sources and observers, and also expanded through 
second order in e. Two main results emerge: i) Higher-order perturbations secularly 
grow at large source-observer distances, and ii) some higher-order metric perturbations 
produce observable effects with angular dependencies which are completely different 
from those associated with first-order effects. 

The first of these statements can be understood directly from (56), which shows 
that the second-order metric perturbation involves the first-order metric perturbation 
hij via a double integral of the “energy density” hij hij. This density is non-negative, 
so its integrals—and therefore the metric itself—grow wherever hij ^0. As a 
consequence, second-order optical effects can be much more important on large scales 
than naive estimates might suggest. 

This can be interpreted as a kind of memory effect when the curvature is 
significant only for a short time (i.e., for gravitational wave bursts). Standard 
assumptions have long been known to imply that at first order, bursts can exhibit 
a displacement-type memory, but no “velocity memory;” initially-comoving particles 
remain comoving after a gravitational wave has passed. We show that this picture 
changes at second order in the gravitational wave amplitude. The secularly-growing 
portion of the second-order metric physically corresponds to a finite kick imparted 
to initially-comoving test particles, a nontrivial velocity memory. Even a small effect 
of this sort can produce significant displacements over sufficiently long times. If a 
particular burst is characterized as N oscillations with angular frequency w and strain 
amplitude e, the associated nonlinear effects on observables are fractionally of order 
e^N{ijjr) over lengthscales of order r 3> Nuj~^. 
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The analogous scaling is different for continuous gravitational waves which 
maintain their amplitudes over all relevant lengthscales. Second-order effects are then 
shown to grow like where N ^ ujr now denotes the number of gravitational 

wave cycles between a source and its observer. This number is typically of order unity 
or less for gravitational waves intended to be observed using standard interferometer 
designs such as LIGO, so nonlinear effects are negligible in those cases. Much larger 
values of N can arise in pulsar timing, however. Consider, for example, a gravitational 
wave with 10127: = 300 nHz and with an approximately constant amplitude between 
the earth and a pulsar where r ^ 10 kpc. Second-order terms in this case are then 
amplified by ~ 10^^. Although this factor is large, multiplying it by a realistic 
strain magnitude results in an O(e^iV^) effect which would still be challenging to 
detect. 

Regardless of the precise type of gravitational wave under consideration, we 
show that the growing higher-order metric perturbations also have nonvanishing 
traces. While the familiar trace-free property of TT gauge can be extended to all 
odd-order metric perturbations, even orders generically acquire finite traces. This 
distinction physically results in different azimuthal dependencies for the optical effects 
associated with, e.g., first and second-order metric perturbations. Focusing only on 
the aforementioned second-order terms which grow at large distances, many second- 
order observables in general relativity appear similar to linear observables, but with 
hij replaced by an “effective metric perturbation” which possesses a small nonzero 
trace. This trace has optical effects similar to those associated with breathing-type 
polarization modes in other theories of gravity. It differs, however, in that nonlinear 
perturbations in general relativity do not represent physically independent degrees of 
freedom; such terms are entirely determined by the ordinary -|- and x polarization 
modes. Furthermore, this effect arises only for nonlocal measurements. Both 
distinctions can be observationally subtle, however, and might complicate efforts [46] 
to constrain alternative theories of gravity via the presence of additional gravitational 
wave polarization modes. 
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